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ABSTRACT 

Saturn is the only known planet to have coorbital satellite systems. In the present work 
we studied the process of mass accretion as a possible mechanism for coorbital satellites 
formation. The system considered is composed of Saturn, a proto-satellite and a cloud of 
planetesimals distributed in the coorbital region around a triangular Lagrangian point. The 
adopted relative mass for the proto-satellite was 10^® of Saturn's mass and for each plan- 
etesimal of the cloud three cases of relative mass were considered, lO^^^, lO^^^ and lO^^^ 
masses of Saturn. In the simulations each cloud of planetesimal was composed of 10'^, 5 x 10'^ 
or 10'' planetesimals. The results of the simulations show the formation of coorbital satellites 
with relative masses of the same order of those found in the saturnian system (10^^'^ - 10~^). 
Most of them present horseshoe type orbits, but a significant part is in tadpole orbit around 
Li or L5. Therefore, the results indicate that this is a plausible mechanism for the formation 
of coorbital satellites. 
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1 INTRODUCTION 

Coorbital systems are those in which at least two bodies share 
the same mean orbit. Coorbital objects that librate around the 
Lagrangian points L4 or Ls are said to be in tadpole orbits, while 
those that librate around L4, L3 and L5 are said to be in horseshoe 
. type orbits. Although Lagrange have described the motion of 
' bodies around these equilibrium points in 1788, when he published 
his Analytical Mechanics, only in 1906 a body showing this kind 
of motion was discovered by Max Wolf in Heidelberg. He found 
an asteroid librating around the point L4 of Jupiter in the system 
Sun- Jupiter. This asteroid is (588) Achiles and at the moment there 
are more than 3200 know asteroids coorbiting with Jupiter (http: 
www.cfa.harvard.edu/iau/lists/JupiterTrojans.html). Apart from 
Jupiter there are other planets with coorbital bodies. In 1991 the 
asteroid 5261 Eureka was discovered librating around the Ls of 
Mars (Innanen, 1991). Later three more asteroids were discovered: 
1998 VF31 (Tabachnik & Evans, 1999), 1999 UJ7 (Connors et al., 
2005) and 2007 NS2 (http: 

www.cfa.harvard.edu/iau/lists/MarsTrojans.html). Neptune has six 
coorbital asteroids. All of them librating around its Lagrangian 
point L4 (Zhou et al., 2008). 

On the other hand, all the coorbital satellite systems known 
are around the planet Saturn (Table l).The satellite Thetys has two 
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coorbital satellites. Telesto around its I/4 point and Calypso around 
its L5 point. They were discovered in 1980 by observations made 
from Earth (Veillet, 1981; Reitsema, 1981). The satellite Dione 
also has two coorbital satellites. Helene around its L4 point and 
Polydeuces around its is point. Helene was discovered in 1980 by 
observations made from Earth (Lecacheux et al., 1980; Reitsema 
et al., 1980) and Polydeuces was discovered through images from 
the Cassini spacecraft (Murray et al., 2005; Porco et al., 2005). In 
1981, images from Voyager not only confirmed the existence of the 
coorbital system Janus and Epimetheus but also determined their 
masses and orbital elements. In a suitable rotating coordinate sys- 
tem, both satellites perform horseshoe orbits and reach a distance 
of only 15000 km from each other (Yoder et. al, 1983). 

The satellites Thetys and Dione are of the same size and have 
mass of the order of 10~® Saturn's mass (Table 2). They also 
have another interesting common feature. They are in mean mo- 
tion resonance with a smaller interior satellite. Mimas-Thetys are 
in 4:2 inclination type resonance while Enceladus-Dione in a 2:1 
eccentricity type resonance. However, Mimas and Enceladus do 
not have any coorbital satellite known until now. A possible ex- 
planation for that was given by Mourao et al. (2006) and Christou 
et al. (2007). Mourao et al. (2006) studied the stability of hypo- 
thetical satellites coorbital to the satellites Mimas and Enceladus. 
Their results showed that these coorbital satellites are in stable or- 
bits. Then, they explored the stability of these hypothetical coor- 
bital satellites putting Mimas and Enceladus at different values of 
semi-major axis, which they could have occupied in the past, along 
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Table 1. Coorbital Satellites of Saturn jGiorgini et aljfT99^ , Porco et al. (2005), Murray 
et al. (2005),Jacobson et al. (2004)) 



satellite 


relative mass*^ 


Orbit 


Density (g / cm^ ) 


Radius {km) 




Polydeuces 


I~ 5 X 10-13 


L5 Dione 




3.25 


25.8 


Helene 


4.48 X 10"" 


L4 Dione 


1.5 


16 ±4 


14.8 


Telesto 


1.25 X 10~" 


L4 Thetys 


1.0 


12 ±3 


1.3 


Calypso 


6.32 X 10-12 


I/5 Thetys 


1.0 


9.5 ± 1.5 


3.6 


Janus 


3.38e-09 


horseshoe 


0.65 


99.3x95.6x75.6 




Epimetheus 


9.67e-10 


horseshoe 


0.63 


69x55x55 





Mass w.r.t. Saturn 
Libration Amplitude 



Table 2. Satel lites of Saturn which have coorbital satellites in 
tadpole orbits jOiorgini et"ai] il996h ) 





Satellites 




satellite 


relative mass_^ 


density (g j cw? ) 


radius {km) 


Dione 


1.92 X 10-°*^ 


1.469 


562.5 


Thetys 


1.08 X 10"°^ 


0.956 


536.3 



" Mass w.r.t. Saturn 

their orbital migration due to tidal their effects. They found that the 
hypothetical coorbital satellites would always be in stable orbits. 
The only exceptions occurred when Mimas and Enceladus were 
placed at semi-major axis where they were in mean motion first- 
order resonance with each other. The hypothetical coorbital satel- 
lites left their tadpole orbits when Mimas and Enceladus were in 
4:5, 5:6 or 6:7 mean motion resonances. Therefore, if there were 
coorbital satellites of Mimas or Enceladus, they may have lost them 
along their orbital migration, when passing through specific mean 
motion resonances. 

Dermott & Murray (1981a,b) presented a theory on the tad- 
pole and horseshoe type orbits for the case of circular and elliptic 
restricted three-body problem. They also described the coorbital 
motion of Janus and Epimetheus through a numerical study com- 
bined with perturbation theory. The study was made for a third body 
of negligible mass. Then, they generalized some of these results for 
the case when the third body disturbs the other two bodies of the 
system. One of the main results found was a relation between the 
coorbital trajectory and its associated zero velocity curve. 

There are several studies on the origin of the trojan asteroids, 
but not too many on the origin of the coorbital satellites. The origin 
of these satellites and Trojans asteroids are believed to be associ- 
ated to one of the following mechanisms (Yoder et al.,1983; Smith 
et al.,1981): (i) capture by gas drag, (ii) congenital formation or 
(iii) rupture of a parent body. 

Chanut et al. (2008) studied the mechanism of capture by drag. 
They considered several values for the relative mass of the sec- 
ondary body, from 10~^ to lO"'^. They explored the orbital evo- 
lution of planetesimals under orbital decay toward the secondary 
body orbit. There are three possible outcomes along such evolu- 



tion: collision with the secondary body, capture in the coorbital re- 
gion or crossing the coorbital region toward the central body. Once 
a planetesimal is captured in a coorbital trajectory, its amplitude of 
oscillation around one of the triangular Lagrangian points shrinks 
until it is fixed over the Lagrangian point. They also showed that the 
location of the Lagrangian points change according to the density 
of the gas. 

With the numerous discoveries of the extrasolar planets 
Beauge et al. (2007) analyzed the formation of hypothetical Earth- 
type coorbital planets in extrasolar systems. The central body was 
taken as equal to one solar mass and the giant planet with Jupiter's 
mass. They considered two scenarios of formation, one gas-free 
scenario and another rich-gas scenario. In both of them there were 
the formation of a single coorbital of the terrestrial type, but due 
to gravitational instability with other bodies, the accretion process 
was not efficient. 

In this paper we study the mechanism of congenital formation 
of coorbital satellites and in some points this mechanism is similiar 
to adopted in Beauge et al. (2007). According to the modem theo- 
ries of planetary formation, in the beginning of the solar system for- 
mation there was a cloud of gas and dust that settled in a disk. Due 
to several factors the gas dissipated along the time. In the first stages 
the particles were of the size of sub-micrometers up to the order of 
centimeters. This disk of gas and dust gradually generated larger 
bodies through physical collisions leading to agglomeration. Suc- 
cessive encounters among them resulted in a large number of meter 
sized aggregates. From these objects originated planetesimals of 
the order of kilometers. Further encounters between planetesimals 
led to the formation of the planet's progenitors of terrestrial size or 
even ten times this size (Armitage, 2007). The formation of plan- 
etary satellite systems are thought to have followed a very similar 
way (Safronov, 1969; Wetherill 1980, 1990; Hayashi et al. 1985; 
Greenberg, 1989). 

In the present work we consider an intermediate stage of Sat- 
urn's satellite system formation. In this stage all the major satellites 
are almost formed (here called proto-satellites), there is almost no 
gas left, but there is still plenty of small planetesimals in the disk. 
Our model of study is characterized by a central body (Saturn), a 
secondary body (proto-satellite) and a cloud of plantesimals. Since 
it is assumed that the proto-satellite has already a significative mass, 
which is several orders of magnitude larger than the planetesimals, 
it defines (together with the planet) the structure of the phase space 
with the five Lagrangian points. Then, we consider a cloud of plan- 
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etesimals in the coorbital region around one of the triangular La- 
grangian points, under the gravitational influence of Saturn and the 
proto-satellite. The collision between planetesimals are considered 
constructive and through this accretion process larger bodies are 
built in the coorbital region. 

This paper has the following structure. In the next section we 
present the methodology and assumptions adopted in our numerical 
simulations. The results of the simulations are presented in section 
3. Finally, in the last section we present some final conmients and 
our conclusions. 



2 METHODOLOGY 

The dynamical system considered is composed by Saturn (the cen- 
tral body), a proto-satellite and a cloud of planetesimals. The mass 
of the proto-satellite was chosen to be 10~^ of Saturn's mass, 
which is close to the masses of Dione and Thetys (Table 2). It is 
placed in a circular orbit with the semi-major axis equal to that of 
Thetys. The mass of the proto-satelite might vary along the simula- 
tion due to some collisions with planetesimals, but its small growth 
do not affect significantly the dynamic of the planetesimals around 
the Lagrangian points (Fleming and Hamilton, 2000). 

The planetesimal clouds were always initially randomly dis- 
tributed in a sector around I/4 or L5. The sector is delimited by 
an arc of 80deg, centered on the Lagrangian point, and the ex- 
treme orbital radii of the largest horseshoe orbit (Figure 1). The 
half width of the largest horseshoe orbit is given by Dermott & 
Murray (1981a) as 

Arhorsc = ^fJ-l^^a2 (1) 

where fi2 and 02 are the relative mass and the semi-major axis of 
the proto-satellite, respectively. Studying the overlap of first order 
mean motion resonances. Wisdom (1980) found an expression that 
gives the width of a chaotic region from the planar circular orbit 
of a secondary body. Therefore, it is important to have in mind 
that just outside the coorbital region (interior and exterior) there is 
a chaotic region. Particles in this region usually tend to have close 
encounters with the proto-satellite and get large eccentricities, leav- 
ing the neighborhood. Since we are interested in the formation of 
the small counterparts of coorbital satellite systems, in general, we 
considered clouds of planetesimals with the total mass of the order 
of the smaller coorbital satellites, given in table I. We also consid- 
ered that all planetesimals have the same initial mass, mp-. We 
performed simulations for different clouds of planetesimals. We 
adopted clouds of 10"', 5 x 10'' and 10* planetesimals each one 
for three different values of mp. , 10~^^, 10"^^ and 10~^* masses 
of Saturn. 

The simulations were made through numerical integrations us- 
ing the Burlish-Stoer integrator from the package Mercury (Cham- 
bers, 1999). The length of the integration, r (in orbital periods 
of the proto-satellite), varied according to the value of mp.. For 
rup. = 10~^^ it was used r — 10^, for mp. = 10^^'^ it was used 
r — 10® and for rrip. — 10~^* it was used t = lO'^. That is neces- 
sary due to the gravitational attraction between the planetesimals, 
which is smaller as smaller are their masses, leading to a longer 
time for the system to evolve. 

For each cloud of planetesimals with a specific value of mp . 
we performed at least two independent simulations. One placing it 
around L4 and the other around L5. We also tested the effect of 
Saturn's oblateness by running extra simulations with the inclusion 
of the J2 term in the gravitational potential of the central body. 




Figure 1. Schematic location of the sectors around L4 and L5 were the 
planetesimals are initially distributed. Each sector is delimited by an arc of 
80 dcg, centered on the Lagrangian point, and the extreme orbital radii of 
the largest horseshoe orbit. The half width of the largest horseshoe orbit, 
Arjiorsc. is given by Equation (1). 

The planetesimals in all simulations were initially placed in circular 
orbits. 



3 NUMERICAL SIMULATIONS 

Along the numerical simulations the planetesimals collide among 
them generating lager planetesimals. The collisions are always con- 
sidered constructive, i.e. inelastic collisions. By the end of the in- 
tegration time there are only a few (at most four) large planetesi- 
mals left. The final masses of the remaining planetesimals are rep- 
resented by TUpj, . 

In order to have an idea of the time evolution of these systems 
we present a sample of snapshots for two representative simula- 
tions. Figure 2 shows the dynamical evolution of a cloud of 10* 
planetesimals with mp. — lO"'^^, initially distributed around L4. 
In each frame are shown the values of the osculating semi-major 
axis and eccentricities of the remaining planetesimals at a given 
moment. The orbital radius of the proto-satellite is considered as 
unity. The color grade indicates the mass of the planetesimal. 

Since we are interested in the formation of small coorbital 
satellites, the best candidates should be the remaining planetesi- 
mals with low eccentricity and semi-major axis inside the coor- 
bital region. So we included in the plots the dashed curves that 
indicate the coorbital region defined by the largest horseshoe width 
(1 ± Arhorsc). In general, the planetesimals that leave the coorbital 
region present chaotic behavior and collide with the proto-satellite. 

From the plots of Figure 2 we notice that the eccentricities of 
the planetesimals increase as they start to grow. However, they are 
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limited to low values (e < 0.007). After about 1068 years there 
are only three planetesimals left and it takes about a thousand years 
to reach a single planetesimal. That last one has a mass of mp^ = 
S.UTOB - 10 

In Figure 3 we present the dynamical evolution of a cloud of 
10'^ planetesimals with rrip^ = 10~^*, initially distributed around 
1/5 . The evolution is very much similar to that presented in Fig- 
ure 2. However, at least two major differences should be pointed 
out. The eccentricities do not get much larger then 0.001 and the 
time scales are much longer then in the previous case. These two 
features are naturally due to the lower values of the initial mass 
of the planetesimals, rrip^, and the total mass of the cloud of the 
planetesimals. Larger masses produce stronger gravitational inter- 
actions accelerating the evolution of the system and increasing the 
values of eccentricities. 

The evolution of all the other numerical simulations we have 
performed are similar to these ones. However, the plots presented in 
Figures 2 and 3 do not guarantee that the remaining planetesimals 
are in coorbital orbits with the proto-satellite. Therefore, we ana- 
lyzed the evolution of the angle 9 between the remaining planetesi- 
mals and the proto-satellite, called libration angle. This angle is de- 
fined as the difference in mean longitudes, 6 = Ap — A2, where \p 
and A2 are the mean longitudes of the planetesimal and secondary 
body respectively. A representative sample of such evolution is pre- 
sented in Figures 4 to 7. In each figure we present two plots. The 
top plot shows the evolution of 9 and the bottom one shows the evo- 
lution the of planetesimal's mass. All plots cover the whole integra- 
tion time and they are shown in a logarithmic scale. The temporal 
evolution of the libration angle confirmed that in all simulations 
the remaining planetesimals are in coorbital trajectories. The plots 
show a complete variety of evolutions. Planetesimals that started 
in tadpole orbits around L4 or L5 and ended in horseshoe orbits 
(Figures 5 and 6). Planetesimals that started in tadpole orbit around 
L4 and ended in tadpole orbit around L5 (Figure 4) or in tadpole 
orbit around L4 (Figure 7). The evolutions are not always smooth. 
There are cases of tadpole orbits that "jump" from L4 to L5 and 
then again back to L4 (Figure 7). And others that "jump" from L5 
to 1/4 and then to a horseshoe orbit (Figure 6). 

The changes in the amplitude of oscillation of the coorbital 
planetesimals along their evolution are associated to the collisions 
with other planetesimals. Several of these changes can clearly be 
seen by comparison between the top and the bottom plots of each 
figure. For example, in Figure 6 the last visible step in the growth of 
the planetesimal coincides with the change from tadpole to horse- 
shoe orbit. 

A summary of the whole set of simulations are presented in 
Tables 3 and 4. We performed 24 simulations neglecting the oblate- 
ness of Saturn (Table 3) and other 18 simulations taking it into ac- 
count (Table 4). In the last two columns the kind of coorbital tra- 
jectory and the final mass of the survivors planetesimals are shown. 
As we can see, in several simulations there were more then one 
planetesimal left. In some cases we followed the simulations for 
longer times. For example, simulation 13a (Table 3) was integrated 
for more then 10^ years and no extra collision occurred along the 
second half of such integration time. In the case of simulation 14a 
(Table 3), the integration was extended to one million years and 
there were still two planetesimals left. One in horseshoe orbit and 
other in tadpole around L5 . 

From Table 3 we have that 2/3 of the survivors are left in 
horseshoe trajectories, while the remaining 1/3 are equally divided 
in tadpole trajectories around L4 and L5. We also note that, in each 
simulation, the total mass of the survivors is very close to the corre- 
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Figure 4. Temporal evolution of the libration angle (top) and the rela- 
tive mass, (bottom) of the remaining planetesimal from the simulation of 
a cloud of 10'* planetesimals with rrip^ = 10~^^, initially distributed 
around L4. After about 60 years the planetesimal reaches its final mass, 
nipj = 5.117 X 10"^". and starts to librate around L5. 




time (years) 

Figure 5. Temporal evolution of the libration angle (top) and the relative 
mass, (bottom) of the remaining planetesimal from the simulation of a cloud 
of 10'' planetesimals with m^^ = 10~^^, initially distributed around L5. 
By the end of the simulation the planetesimal reaches its final mass, rrip^ = 
7.35 X IQ-iO and stays in a horseshoe orbit. 




Figure 6. Temporal evolution of the libration angle (top) and the relative 
mass, (bottom) of the remaining planetesimal from the simulation of a cloud 
of 10'' planetesimals with rrip^ = 10~^^, initially distributed around L4. 
By the end of the simulation the planetesimal reaches its final mass, rrip^ = 
10~*" and stays in a horseshoe orbit. 



COORBITAL SATELLITES OF SATURN: CONGENITAL FORMATION 5 



0.004 r— 

0.0035 - 

0.003 - 

0.0025 - 

0.002 - 

0.0015 - 

0.001 - 

0.0005 - 

— 



0.98 0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 
Semimajor axis 



1.368925 years 



0.004 
0.0035 

0.003 
0.0025 

0.002 
0.0015 

0.001 
0.0005 



0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 
Semimajor axis 



16.42710 years 



0.004 
0.0035 



0.002 
0.0015 



0.98 0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 
Semimajor axis 



34.22313 years 



0.004 
0.0035 

0.003 
0.0025 

0.002 
0.0015 

0.001 
0.0005 



0.98 0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 
Semimajor axis 



65.70841 years 



0.004 r 

0.0035 - 

0.003 - 

0.0025 - 

0.002 - 

0.0015 - 

0.001 - 

0.0005 - 

— 



0.98 0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 
Semimajor axis 



0.004 
0.0035 

0.003 
0.0025 

0.002 
0.0015 

0.001 
0.0005 




• 



0.98 0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 
Semimajor axis 



1003.422 years 



0.004 r 

0.0035 - 

0.003 - 

% 0.0025 - 

I ^-^^^ ' 

^ 0.0015 - 

0.001 - 

0.0005 - 

— 



0.98 0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 
Semimajor axis 



1006.160 years 



0.004 
0.0035 
0.003 
^ 0.0025 

I ^-^^^ 
^ 0.0015 
0.001 
0.0005 




0.98 0.985 0.99 0.995 1 1.005 1.01 1.015 1.02 
Semimajor axis 



Figure 2. Dynamical evolution of a cloud of 10* planetesimals with mp^ = 10"^"^, initially distributed around L4. In each frame are shown the 
values of the osculating semi-major axis and eccentricities of the remaining planetesimals at a given moment. The color grade indicates the mass of the 
planetesimal. The dashed curves indicate the coorbital region defined by the largest horseshoe width (1 ± Arijorso)- 



spending total mass of initial cloud of planetesimals. The difference 
between these two values is the mass of planetesimals that collided 
with the proto-satellite. 

Comparing the simulations presented in Table 3 with those 
presented in Table 4, we can note that the oblateness of Saturn did 
not introduce any significant change in the results as a whole. 

Beauge et al. (2007) found that the final mass of the formed 
trojans the kind of terrestrial planets did not depend significantly 
neither on the initial mass, nor on the number of the initial pop- 
ulation of planetesimals. In our simulations we do realize that the 
number of planetesimal and its initial mass play an important role 
in the final outcome. 



In Table 3, we can see that the initial populations with small 
mass (10~^^,10~^*) and small number of planetesimal (1000) have 
a larger tendency to yield more than one final body. This effect is 
caused by the weak gravitational attraction when the planetesimals 
are smaller and more distant from each other, consequently few col- 
lisions occur. However, when we have a lot of planetesimals, even 
with small mass, they are closer, then the gravitational attraction is 
larger and collisions become more frequent. 

Although in most cases the final total mass is almost equal to 
the initial value, in one case there is a significant mass loss (Ta- 
ble 3, simulation 3). This occurred due to a collision of a larger 
planetesimal with the secondary body. However, for a similar ini- 
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Figure 3. Dynamical evolution of a cloud of 10"^ planetesimals with rrip^ = 10"^*, initially distributed around L5. In each frame are shown the 
values of the osculating semi-major axis and eccentricities of the remaining planetesimals at a given moment. The color grade indicates the mass of the 
planetesimal. The dashed curves indicate the coorbital region defined by the largest horseshoe width (1 ± Arjujrs;,)- 



tial planetesimal population distributed around L5 there is no mass 
loss, this loss mass might be one statistical fluke, but this has to be 
checked with further studies. 

The temporal evolution of the total number of the remaining 
planetesimals, A'^, in the system gives an idea of the different speeds 
at each stage of the system evolution. In Figure 8 we present a 
representative sample of the temporal evolution of the total num- 
ber of the remaining planetesimals of our simulations. In all the 
cases studied, the evolution of the system can be well represented 
by three stages. In the first stage there is a "cold" cloud, the plan- 
etesimals are in near circular orbits and do not collide frequently. 
Then, along the time the gravitational interaction among the plan- 



etesimals "heat" the cloud, the trajectories get some eccentricity 
and the rate of collisions increase, producing larger planetesimals. 
Larger planetesimals excite the cloud even more and also increase 
their sphere of influence, resulting in larger cross sections and faster 
growing. When most of the planetesimals collided, there are a few 
left, then the collision rate decreases significantly and the evolution 
of the system slows down. 

A comparison of the curves presented in Figure 8 shows that 
the rate of decay of the total number of planetesimals is strongly 
dependent on the total mass of the initial cloud of planetesimals. 

In order to infer possible effects due to mean motion reso- 
nances with nearby satellites on the mechanism of coorbitals for- 
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Table 3. Set of simulations without Saturn's oblateness 



Simulation 


Number of 
plantesimals 


Location of 
initial distribution 


Integration 
time (years) 




Trajectories of 
the survivors'* 


final masses 

nip J, 


1 


10000 


L4 


10005 


10-14 


H 


10-1" 


2 


10000 


L5 


20238 


10-" 


H 


10-1" 


3 


10000 


U 


1068 


10-13 


Tl, 


5.117 X 10-1'' 


4 


10000 


L5 


99 


10-13 


H 


9.999 X 10-1° 


5 


10000 


U 


51 


10-12 


H 


9.969 X 10-^ 


6 


10000 


is 


50 


10-12 


H 


9.974 X 10-^ 


7 


5000 


Li 


51000 


10-14 




2.297 X 10-11, 2.703 x IQ-" 


8 


5000 


Ls 


51000 


10-" 




2.53 X 10-11, 2.47 x 10-" 


9 


5000 


Li 


3325 


10-13 


H 


4.961 X IQ-i" 


10 


5000 


L5 


989 


10-13 


H 


4.943 X 10-1" 


11 


5000 


La 


52 


10-12 


H 


3.602 X 10-^ 


12 


5000 


L5 


61 


10-12 


Tl, 


3.369 X 10-3 


13a 


1000 


Li 


51000 


10-" 


H, H, Tl^ 


3.45 X 10-12, 1.94 X 10-12, 4.61 x IO-12 


13b 


1000 


Li 


51000 


10-" 




2.54 x 10-12, 2.21 X 10-12, 5.25 x IO-12 


14a 


1000 


L5 


51000 


10-" 




3.35 X 10-12, 9.6 X 10-13, 4.92 x IO-12, 7.7 x IQ-" 


14b 


1000 


L5 


51000 


10-" 


H, H, H, Tl^ 


1.51 X 10-12, 5.03 x 10-12, 2.14 x IO-12, 1.33 x IO-12 


15a 


1000 


Li 


5100 


10-13 




5.38 X 10-11,4.62 x 10-" 


15b 


1000 


Li 


5100 


10-13 


T'l4,H 


4.24 X 10-11, 5.76 x 10-" 


16a 


1000 


L5 


5100 


10-13 




1.7 X 10-12,4.32 X 10-11,5.51 x IQ-" 


16b 


1000 


L5 


5100 


10-13 




4.28 X 10-11, 5.72 x 10-" 


17a 


1000 


Li 


510 


10-12 




6.99 X 10-1", 2.52 x lO-"' 


17b 


1000 


Li 


229 


10-12 


H 


9.32 X 10-1° 


18a 


1000 


L5 


192 


10-12 


H 


7.35 X 10-1° 


18b 


1000 


L5 


277 


10-12 


H 


10-3 



Horseshoe (H), Tadpole around Li (T^^) and Tadpole around L5 (Tlj) 

mation, we performed simulations considering Mimas and Thetys 
in 4:2 resonance of the type-inclination. We considered a sys- 
tem composed of an oblate central body (Saturn), two satellites 
(Mimas-Thetys) which are in mean-motion resonance 4:2 and a 
cloud of planetesimals. The simulations were made for six cloud 
of 1000 planetesimal considering in each case rUp- equal to IO-12, 
10-13 or 10-14 around L4 and L5. It is important to have in mind 
that the previous simulations were in a planar system, and when 
we considered Mimas and Thetys in a 4:2 inclination-type reso- 
nance, the problem becomes three-dimensional. In general, these 
new simulations did not show any significant difference in com- 
parison with the simulations presented in table 4. In Figure 9 we 
present a comparison of the temporal decaying number of the plan- 
etesimals for the simulations of a cloud of 103 planetesimals ini- 
tially located around L4 and with nip. = IO-12. It shows that only 
in the final stages (after reducing the number of planetesimals to 



about 30%) the results present some notable difference. The results 
under resonant perturbation evolves more slowly than the system 
without such perturbation. Probably, this is only due to the fact that 
the system under resonant perturbation is in a three-dimensional 
space, which reduces the chances of encounters between the plan- 
etesimals in comparison with the planar system case. 



4 FINAL COMMENTS 

In the present work we present results of numerical simulations in 
order to study the viability of the congenital formation approach 
as a mechanism for the origin of the coorbital satellites of Thetys 
and Dione. We considered that the collisions are always inelastic 
and make the growth of the planetesimals. The results show that 
almost all the initial mass of the clouds are converted in one or a 
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Table 4. Set of simulations with Saturn's oblateness 



Simulation 


Number of 
plaiitesimals 


Location of 
initial distribution 


Integration 
time (years) 




Trajectories of 
the survivors'* 


final masses 


1 


10000 


Li 


51000 


10-" 


H, H 


1, 612 X 10-11, 8.388 x IQ-" 


2 


10000 


L5 


14672 


10-" 


H 


10-1" 


3 


10000 


U 


654 


10-13 


H 


8.666 X 10-1'^ 


4 


10000 


L5 


218 


10-13 


H 


9.822 X 10-1° 


5 


10000 


U 


55 


10-12 


H 


9.382 X IQ-'J 


6 


10000 


is 


59 


10-12 


H 


9.923 X 10-'-* 


7 


5000 


Li 


51000 


10-" 




2.646 X 10-11, 2.354 x IQ-" 


8 


5000 


Ls 


51000 


10-" 




2.119 x 10-11, 2.881 X 10-11 


9 


5000 


Li 


2308 


10-13 


H 


4.926 X 10-1" 


10 


5000 


L5 


1663 


10-13 


H 


5 X 10-"' 


11 


5000 


La 


95 


10-12 


H 


4.445 X IQ-'J 


12 


5000 


Ls 


49 


10-12 


H 


4.333 X lO-'J 


13 


1000 


Li 


51000 


10-" 


H, H.H, Tl, 


3.18 X 10-12, 3.45 X 10-12, 1.79 x 10-12,1.58 x IO-12 


14 


1000 


Ls 


51000 


10-" 


H, n,TL, 


4.39 X 10-12, 2.12 X 10-11, 3.49 x IQ-" 


15 


1000 


La 


29768 


10-13 


H 


6.78 X 10-11 


16 


1000 


Ls 


51000 


10-13 




4.37 X 10-11, 5.63 x 10-" 


17 


1000 


La 


218 


10-12 


H 


9.86 x 10-1" 


18 


1000 


Ls 


173 


10-12 


H 


9.99 x 10-1" 



Horseshoe (H), Tadpole around La (Tl^) and Tadpole around L5 (Tlj) 
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Figure 7. Temporal evolution of the libration angle (top) and the relative 
mass, (bottom) of the remaining planetesimal from the simulation of a cloud 
of 5 X 103 planetesimals with mp^ = 10" 1^, initially distributed around 
L4. After about 51 thousand years the planetesimal reaches its final mass, 
mp^ = 2.703 X 10—11 and starts to librate in a tadpole orbit around L4. 

few planetesimals in coorbital motion with the proto-satelhte. The 
final masses are of the same order of the small coorbital satellites 
Helene, Polydeuces, Telesto and Calypso. Therefore, we found that 
this is a promising mechanism for this purpose. 



The distribution of librational amplitudes of the final bodies 
that remain in tadpole orbits varies from 4 up to 75 degrees. The 
real coorbitals have a libration amplitude of which varies from 1.3 
up to 25.9 degrees. However, some mechanisms could reduce the 
amplitude of libration, for example, the gas drag (Chanut et al., 
2008), mass accretion or radial migration of the secondary body 
(Fleming and Hamilton, 2000). In this paper we do not consider 
these mechanisms, but might be included in future works. 

The effects of mean motion resonances of nearby satellites 
shows to be of no significative relevance for the present model of 
congenital formation of coorbitals. 

The results presented here might look like obvious. However, 
in Beauge et al (2007) it was found that a growing planetary embryo 
generates a chaotic region around the Lagrangian point, inhibiting 
additional accretion and setting a limit in the final mass of the coor- 
bital body. They studied the possibility of congenital formation of 
Earth like planets in coorbital orbits of exoplanetary systems. In 
their work they adopted systems like the Jupiter-Sun system, with 
relative mass of 10-3, and did not manage to form coorbital plan- 
ets with mass larger than 0.6iV/® This phenomena does not occur in 
the present work because the results of the simulation are strongly 
dependent on the relative mass of the secondary body (Izidoro et 
al., 2010). 

It is important to comment that in the present congenital model 
the final coorbitals formed are of similar masses. However, Dione 
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Figure 8. Temporal evolution of the total number of remaining planetesi- 
mals. It is shown the data from simulations of clouds of 10'^ (top), 5 X 10'^ 
(middle) and 10^ (bottom) planetesimals around L4 with initial masses 



10" 



10" 



' and rUr, 



10" 



has one coorbital, Helene whose mass is about two orders of mag- 
nitude higher than the other coorbital, Polydeuces. Therefore, the 
model as considered in the present work is not able to explain by 
itself alone such configuration. Consequently, there is still plenty 
of work to be done. 
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